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The conventional color-singlet model was challenged by the recent data on quarkonium 
production. Discrepancies in production rates were observed at the Tevatron, at LEP, 
and in fixed-target experiments. The newly advocated color-octet mechanism provides 
a plausible solution to the anomalous quarkonium production observed at the Tevatron. 
The color-octet mechanism should also affect other quarkonium production channels. In 
this paper we will summarize the studies of quarkonium production in Z° and T decays. 

1. Introduction 

Charmonium was first discovered more than twenty years ago, which was be- 
lieved to be a pair of charm and anti-charm quarks bounded by strong QCD forces. 
The production and decay of quarkonia had then become very interesting subjects 
in QCD. The production of quarkonia was first formulated by the color evaporation 
modelEl. It states that the production rate of a particular quarkonium state H is 
a certain fraction fu of the open charm pair production with invariant mass m c g 
between 2m c and 2mr), i.e., for example, 

f 2m ° ct(cc) 

OW = U X / -j dm cS . (1) 

The model has the predictive power that /#'s are process independent, and can 
reproduce the y/s dependence of the total cross section. However, it is unable to 
predict the relative rates for different spin-orbital states. 

Tha color evaporation model was later superseded by the color-singlet model 
(CSM)0. This model is based on the assumption that the spin-orbital angular 
momentum and the color of the asymptotic physical state is the same as the point- 
like cc pair produced in the short distance region. Thus, the amplitude involving 
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a physical charmonium state is proportional to the amplitude involving the point- 
like cc pair in the same angular momentum state with cc in a color-singlet state 
and having small relative momentum. The proportionality constant is given by the 
wavefunctions or derivatives of wavefunctions evaluated at the origin, which are 
obtained by solving the Schroindger equation of the c and c system. 

The CSM is simple with high predictive power. However, it often gives quarko- 
nium production rates below those measured from the fixed-target experiments. 
The contrast becomes much prominent when measured rates of prompt tp and tp' in 
CDF experiment^ exceed those predicted by the CSM by orders of magnitude at the 
large transverse momentum region. The inadequacy of CSM led to two important 
theoretical developments in quarkonium production^. The first one is the idea of 
gluonEl and heavy quark fragmentationB. Though these fragmentation processes are 
higher order corrections in a s , they have been shownQ to be more important than 
the lowest order QCD production^ in the large pr region. The dominant contri- 
bution to tp production would be gluon fragmentation g* — > XcJ followed by the 
radiative decays XcJ — > "07- For tp' the dominant contributions come from gluon 
and charm quark fragmentation. However, ihe tp data showed that the tp fed down 
from XcJ decay is not the major source of tpa. Moreover, the tp' data was still orders 
of magnitude above the prediction even including fragmentation contributions. 

The second important development is the color-octet mechanismH, which in- 
cludes higher order corrections in v 2 , where v is the relative velocity of c and c 
inside the charmonium. Within Non-Relativistic QCD (NRQCD)li3 formalism of 
the quarkonium, the physical state, e.g. tp, is not solely a cc pair in 3 Si color- 
singlet state, instead, a superposition of all Fock states: 

|0) = \cc( 3 S 1 ,l))+0(v)\cci 3 P J ,8)g)+0(v 2 )\cc( 3 S 1 ,8orl)gg) 

+ O(v 2 )\cc( 1 S ,8)g)+O(v 2 ) \cc( 3 D,,,, 8 or l)gg) + 0(v 3 ) , (2) 

where the angular momenta of the cc pair in each Fock state is labeled by 2S+1 Lj 
with a color configuration of either 8 or L Similar expressions exist for XcJ states. 
Production of tp can now be via higher Fock states with emission or absorption 
of dynamical gluons. The NRQCD factorization formula for the production of a 
quarkonium state H is given by 

da(H + X) = ^2 da(cc(n) + X){Q%) , (3) 

n 

where da is the short distance factor producing the cc pair in an intermediate state 
n (n is a collective notation for angular momentum and color), and (O^) is the 
matrix element describing the transition of the state n into the physical state H. 
The relative importance of each term in (||) depends on the order of a s in da 



and v 2 in (£>%}. The 
The best known example 



essentially the leading term in the expansion (J3J). 
I of color-octet mechanism is the gluon fragmentation 
into a cc( 3 Si,8) intermediate state, which then evolves nonperturbatively into the 
physical tp or tp' . This is so far the best plausible solution to the anomalous tp and 
tp' production at the large px region in hadronic collisions. 



Quarkonium Production at Z° and in T(1S) Decay 3 



If the color-octet mechanism is correct, it would have significant impact on 
other pr odu ctio n modes, e.g photoproductiorilj. hadrmproductionO, production 
in Z°ESE3, TtH, and B decayai3, and at e + e~ collisionaLj. While othermodes are 
summarized by other speakers, we concentrate on the production in Z°Ej and TllJ 
decays. 



2. ijj Production in Z° Decav 

Though the major source of ip in Z° decay comes from B meson decays, the LEP 
detectors are able to separate the prompt ip from the ip coming from B decay. The 
recent data on prompt ip production showed a factor of almost 10 larger than the 
prediction by the CSM. Therefore, it would be interesting to examine the effects of 
the color-octet mechanism in this production mode. 

Let us first summarize a few important color-singlet_processes. The leading order 
color-singlet processes are Z — > ■059" an d Z ~ * ipcc£^. Although both processes 
are of order a 2 s , the latter is two orders of magnitude larger. The latter process 
has been interpreted as a fragmentation processcl in which the cc pair was first 
produced almost on-iihcll from the Z decay and then followed by the fragmentation 
of c or c into the ipu. Therefore, the process Z — > ipec is not suppressed by the 
quark propagator effect of order 1/Mz as it does in the process Z — > ipgg. In the 
fragmentation approximation, the energy distribution of ip in the process Z — > ipec 
is given byu 

(Z -> iP{z) cc) « 2 T(Z -> cc) x D c ^(z) , (4) 



dz 



where 



D ^ iz) = 1 ^&0<!M^r (16 _ 32s+72 ^_ 32 ,3 +fc4) , (5) 

and z — 2E^,/Mz with denotes the energy of ip. Both processes involve the 
same matrix element (Of ( 3 Si)) whose value can be extracted from the leptonic 
width to be about 0.73 GeV 3 . Numerically, the widths for Z — > ipgg and Z — * ipec 
are about 6 x 10~ 7 GeV and 7 x 10 -5 GeV, respectively. There is also ahigher order 
color-singlet process Z — > qqg* followed by the gluon fragmentatior£3 g* —> ipgg. 
The branching ratio was estimated to be of order 10~ 6 onlySii. However, the recent 
results from LEPcJ showed a branching ratio of order 10 -4 for prompt ip production, 
which is well above all the predictions from the color-singlet model. 

The leading order color-octet process is of order a s given by the process Z — * tpg, 
for which one of the Feynman diagrams is shown in Fig. 0(a). Unfortunately, the 
short-distance factors are suppressed by at least one power of 1/Mz- Numerically, 
the width for Z — > ipg is of order 10 -7 GeV which renders this process useless at 
the Z resonance. 

The dominant color-octet process actually begins at the order a 2 s in the process 
Z — > ipqq, for which one of the Feynman diagrams is shown in Fig. 0(b). Here q 
represents u, d, s, c, or b. The energy distribution of ip for the process Z — > ipqq is 
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Fig. 1. Representative Feynman diagrams for (a)Z — > j/jg, and (b) Z — > tpqq. 



calculated, in the limit m„ = 0, to befe? 



dT 

— {Z -> ip(z)qq) 



a " {2mc) -T(Z^ q q)^ {3Sl)) 
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where z = 2E^/M Z , £ = Ml/M%, and z L = (z 2 - 4^) 1 / 2 . The physical range of z 
is 2y/£ < z < 1 + ^Integrating Eq. (|^) over the physical range of z gives the decay 
width of Z 



T(Z ^ j;qq) _ aj(2m c ) (Og^)) 
T(Z -> to) 



2 Li, 
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5(l-e 2 )-2Cloge+(l + 2 x (7) 
21og(l + 01oge + 31oge + log 2 e 



where 1^2(3;) = — J Q * log(l — i) is the Spence function. Numerically, the above 
ratio is 2.2 x 10 -4 . Summing over all the quark flavors (q — u, d 1 s, c, 6), we obtain the 
decay width T(Z — > ^><7f) ~ 3.8 x 10~ 4 GeV and the branching ratio Br(Z — > 
■09?) — 1.5 x 10~ 4 . Assuming the dominant prompt ip production process to be 
Z — > (75(7* with — > ^ + X, in which according to color-singlet model the off- 
shell gluon fragments into a ip plus two perturbative gluons, DELPHlEa obtained 
the limit Br(Z — > j* ^ + 1) < 4.1 x 10~ 4 . Since the event topology of 

our color-octet process is similar to this one, this limit should also be valid for our 
color-octet process. Therefore, our result is consistent with this data. 

In Fig. we compare the energy distributions of ip coming from the most impor- 
tant color-octet process Z —* ipqq using Eq. (Jfy and the most dominant color-singlet 
process Z — > ipcc using Eq. (§. The comparison in Fig. |^ shows a very spectacular 
difference between the color-octet and color-singlet processes. The spectrum for the 
color-octet process is very soft with a pronounced peak at the lower z end, while 
the spectrum for the color-singlet process is rather hard because of the nature of 
heavy quark fragmentations. 

Another possible test to distinguish the charm quark fragmentation Z — > ipcc 
and the color-octet process Z — > tpqq is using the polarization of ip. The polarization 
of ip can be measured by the angular distribution of the lepton pair in the rest 
frame of ip. The charm quark fragmentation process predicts a polarization ratio 
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Fig. 2. Energy spectrum dT/dz of tfi from the color-octet process Z — > ?/j<jr<j and the color-singlct 
process Z — > i/icc. 



L/(L + T) « 31%, while the dominant color-octet process gives a ratio L/(L + T) » 
20%. Nevertheless, the small production rate of prompt r/> renders this test rather 
difficult . 

Concluding this section the color-octet process Z — ► contributes dominantly 
to the prompt ^ production, and the resulting branching ratio is consistent with 
the data. In addition, the energy spectrum of tp serves as a useful test to distinguish 
between the color-octet process Z — ► ipqq and the color-singlet process Z — * ^cc. 

3. "0 Production in T DecayEI 

The available experimental data on charmonium production in T decay are listed 
as follows, 



Br(T-> ip + X) 



= (1.1 ± 0.4) x 10~ 3 CLEO 24 , , 

< 1.7 x 10~ 3 CrystalpBall 

< 0.68 x 10~ 3 ARGUSEa. 



(8) 



Such large branching ratios by CLEO can hardly be explained by the CSM because 
the lowest order color-singlet processes are of order a® and there is no quantitative 
calculation of such existing in the literature. Only until a few years ago an indirect 
tp production via the intermediate physical \c states was performed by TrottierEZI. 
This process is of order a\, i.e., l/a s larger than the color-singlet processes. Never- 
theless, the branching ratio obtained is still one order of magnitude below the data. 
Therefore, it would be interesting to identify important color-octet processes. 

The NRQCD factorization formulalij for the inclusive charmonium production 
in bottomonium decay, such as T — > ip + X, is given by 



dT(T -► ^ + X) = d^r, 



(T\O m \T)(Oi) 



(9) 
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where dT mn are the short-distance factors for a bb pair in the state m to decay 
into a cc pair in the state n plus anything, where m, n denote collectively the 
color, total spin, and orbital angular-momentum of the heavy quark pair. dT mn 
can be calculated in perturbation theory as a series expansion in a s (m c ) and/or 
a s (mfc). The nonperturbative factor (T|O m |T) is proportional to the probability 
for the bb pair to be in the state m inside the physical bound state T, while (0%) 
is proportional to the probability for a point-like cc pair in the state n to form 
the bound state ip. The relative importance of the various terms in the above 
double factorization formula (^|) can be determined by the order of or v c in the 
NRQCD matrix elements and the order of a s in the short-distance factors dT mn . 
In fact, the calculation by TrottierEZl is an example of this factorization formula. 
An infrared divergence occurs in the leading order calculation of the short-distance 
factor, therefore, in addition to the color-singlet matrix element (0* cj ( 3 Pj)), one 
also needs to include the color-octet matrix element (0 8 cj ( 3 Si)) to absorb the 
infrared divergence. In this case, the introduction of the color-octet matrix element 
is required by perturbative consistency. 

For color-octet processes we first consider the produced cc pair in the color-octet 
1 Sq, 3 Si, or 3 Pj configuration, which subsequently evolves into physical ip described 
by the matrix elements (Of^So)), <0^( 3 5i)), or (0^( 3 Pj)), receptively. The 
matrix element (0% ( 3 Si)) has been extracted from the CDF dataHliM2l, while two 
different combinations of the other color-octet matrix elements have been extracted 
from the CDF datalia and from the photoproduction data by Amundson et a£3. 
Though of much smaller effects, we also consider the contributions by the higher 
Fock state of the color-octet bb pair inside the T associated with the matrix element 
(T|0 8 ( 3 Si)|T), whose value has not yet been determined. An order of magnitude 
of this matrix element can, in principle, be estimated by considering the ratio 

(TjOgfigOrY) (Of( 3 Si)) (vl\ 2 

(TIO^SOIT) (0t( 3 S 1 )) ~ WJ ' 1 1 

which implies (T|0 8 ( 3 Si)|T) w 3 x 1(T 3 GeV 3 . The ratio in @) tells us that 
processes associated with a color-octet cc pair inside the produced ip an d a color- 
singlet bb pair inside the decaying T are much more important than those with a 
color-octet bb pair inside the T and a color-singlet cc pair inside the ip. However, 
such a large value for this matrix element would substantially increase the hadronic 
width of T, which would diminish the leptonic branching ratio to an unacceptable 
level. In order not to spoil the experimental value for the leptonic branching ratio 
and the total hadronic width of T, we necessarily obtain the following bound on 
(T|O s ( 3 ^)|T>: 

(T|0 8 ( 3 ^)|T) « (l.9 +Jg ) x 1CT 4 GeV 3 . (11) 

Therefore, the color-octet processes associated with (T|0 8 ( 3 Si)|T) have very small 
branching ratios, and thus negligible. 



Quarkonium Production at Z° and in T(15) Decay 7 




^ cc ( Lj,8) 



bb ( J S r 1) 



Fig. 3. Some of the contributing Feynman diagrams for the short-distance processes: (a) 
b5( 3 5i, 1) -> 7* -* cc( 2S+1 £j,8) g ; (b) 66( 3 Si,l) ~» 9*99 ~» cc( 3 Si,8)gg; 

The first color-octet process we consider is the one with the short distance 
process bb( 3 Si, 1) — » cc( 2S+1 Lj, 8) +5, which is of order a 2 a s - One of the Feynman 
diagrams is shown in Fig. |3](a). _Using the heavy quark spin symmetry relation 
(Ot( 3 Pj)) ~ (2J+ 1)(C^( 3 P O ))0, the total width from these processes can be 
simplified as 

r„(T^ + X) = (T|0,(35 1 )|T) ( (0?( , So))(1 _ (12) 

3 mtm c ( 



3m 2 



(1-30 6(1+0 2(1 + 3^ + 6^) 



UsingEl (Of^So)) w (O^( 3 P ))/m2 w 10~ 2 GeV 3 , the contribution from the above 
color-octet processes to the inclusive branching ratio BR(T — > V + -X") is oni y 
1.6 x 10 -5 , which is almost two orders of magnitude below the CLEO data (|J). 

The next process we consider is bb( 3 Si, 1) — > ggg* — > cc( 3 Si, 8) + gg, which is of 
order a*. Fig. |3|(b) shows one of the six Feynman diagrams for the bb( 3 S\, 1) pair 
annihilating into three gluons with one of the gluaas converting into the cc( 3 Si,8) 
pair. The partial width is this process is given bycZI 

r b (T^ + *) = 2 {r\o 1 ( 3 s 1 )\r){ot( 3 s 1 )) x (0.571) . (13) 

486mgm^ 

We note_jthat the color-octet piece for the process T — > XcJ + X considered by 
TrottieiO can be obtained from the above Eq. ([l3]) by simply replacing the matrix 
element (Of ( 3 Si)) with (0^ aJ ( 3 Si)). The prediction of the partial width is sensitive 
to the values of the two NRQCD matrix elements, the running coupling constant 
a s , and the heavy quark masses. For convenience we can normalize this partial 
width to the three-gluon width T(T —> ggg) 

r fe (T -> f + X) = na s (OjfSJ) 0.571 
r(T -> ggg) 8 m 3 n 2 - 9 ' 1 ' 
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With (C$( 3 Si)) = 01)14 GeV 3 Bll3, and assuming BR(T -► w BR(T -> 

light hadrons) = 0.92E3, we obtain 

r b (T^y. + x) r b (T^^ + x) 

IWT) = f(T^) BR(T — > ggg) 

w 2.5 x 10~ 4 . (15) 

This prediction is smaller than the CLEO data by merely a factor of 4, and is 
consistent with the bounds from Crystal Ball and ARGUS. 

The above color-octet process can be extended to ip' as well, simply by replacing 
the matrix element (Of( 3 Si)) with (Of I 3 r St)A, whose_yalue has also been extracted 
from the CDF data to be 0.0042 GeV 3 l'Oti. WithS BR(0' -> ip + X) w 57%, 
we obtain a branching ratio of 4.3 x 10~ 5 for the inclusive ip production in the T 
decay that comes indirectly from ip' . One can also consider the processes T — > "fgg* 
followed by g* — > ip (-0') via the color-octet mechanism, and T — > ggj* followed 
by 7* — > ip (ip') in the color-singlet model. Their partial widths are suppressed 
by factors of 8a/(15a.) ~ 0.02 and 32a 2 (Of ( 3 Si)}/(45a! 2 (Of ( 3 Si))) ~ 0.06, re- 
spectively, compared with the width of Eq. ([l3]) . Thus they contribute a branching 
fraction about 2 x 10^ 5 in the inclusive decay T — > ip + X. The indirect contribution 
from the ip' from these two processes is about 6 x 10~ 6 . 

Another process we consider is that the bb pair annihilates into a qq pair via 
the s-channel photon 7* , and then a bremsstrahlung virtual gluon emitted from the 
light quark line splits into an octet cc( 3 Si, 8), which then turns into ip. This process 
is of order a 2 a 2 , which is similar to the dominant color-octet process in the prompt 
ip production in Z° decayEl. The partial width is given by 

r le (T-7*-fl#) _ «2(2m c ) (0*(*Si)) 



T(T -> 7* -> qq) 36 



5(l-n-2^1og^ 



2 u - ( YTt) ~ 2Li2 (tti) _ 2l0g(1 + e) loge 
(i + 2 l. ( l6 ) 



+31og£ + log 2 £ 

with £ = (rric/mb) 2 . Numerically, the ratio in ( pif ) is only 8.3 x 10~ 6 , and so this 
process can be safely ignored. 

Among all contributions, the largest comes from the process shown in Fig. [|(b), 
which has a branching ratio about 2.5_x 10~ 4 . The next largest contribution is the 
indirect process considered by TrottiercZl, which gives a branching ratio of 5.7x 10~ 5 . 
Other comparable processes are (1) the process shown in Fig. ||(a) which has a 
branching ratio of 1.6 x 10~ 5 , (2) the processes T — > ggj* — > ip + X and T — > 
jgg* — > ip + X via the color-octet mechanism which have a combined branching 
ratio about 2 x 10~ 5 , and finally, (3) the indirect contribution from ip' having a 
branching ratio about 5 x 10~ 5 . Therefore, adding up the contributions from all 
these processes, we obtain a branching ratio BR(T — > ip + X) w4x 10~ 4 , which is 
within 2er of the CLEO data. 
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4. Conclusions 

The color-singlet model, which has been popularly used for more than 15 years, is 
inadequate to describe the present data at the Tevatron, LEP, T decay, and fixed 
target experiments. The NRQCD description of quarkonium production provides a 
consistent theoretical framework and implies the so-called color-octet mechanism. 
It can explain the large px ip and ip' data at the Tevatron by fitting a small set of 
parameters of matrix elements. If the color-octet mechanism is correct, it would 
have a strong impact on other production modes. In this talk, we summarized the 
ip production in Z° and T decays. The color-octet mechanism turns out to be very 
important in these modes. The color-octet process Z — > ipqq can explain the data 
on ip production in Z° decay. A more decisive test would be a measurement of the 
energy profile of the ip where a very soft spectrum is predicted by the color-octet 
mechanism. Furthermore, in T decay the color-octet mechanism makes possible a 
new short distance process 6( 3 Si,l) — > ggg* — > ggcc( 3 S\, 8), which can bring the 
theoretical prediction to be within 2a of the central value of the CLEO data. 
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